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l  ABSTRACT 

*  .1  ! 

,  '  :  i  : 

!  , 

The  Krook  kinetic' equation  is  used'as  the  governing  equation 

».  I  .  ’ 

in  ’an  .investigation  of  a  class  of  exterior  flow  problems  involving 
\  *  « 

i  1  t 

flow  past  an  elliptic  cylinder.  A  numerical  method  for  solving  the 

■  «  '  J  >  , 

integral  equation  form  of  the  Krook  equation  for  this  class  of  problems 

i  '  ■  •  ‘ 

is  described,  and  several  sets  of  solutions  are  presented. 
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1 .  INTRODUCTION 


The  Krook  kinetic  equston  [6]  provides  a  useful  characterization 
of  many  flow  problems  in  kinetic  theory.  By  integrating  this  equation 
along  its  characteristics.-,  a  coupled  set  of  nonlinear  integral  equations 
for  density,  velocity  an?  temperature  can  he  obtained.  This  paper 
describes  a  numerical  p- vedure  which  can  be  used  to  solve  these  integral 
equations  for  a  class  exterior  flow  problems  involving  flow  past  an 
elliptic  cylinder. 

A  si nulation  technique  has  been  applied  by  Vogenitz  et.al.  to 
several  exterior  flow  problems,  including  the  leading  edge  problem 
[15]  and  flow  past  cylinders  and  spheres  [14].  However,  the  only  finite 
Knudsen  number  solutions  of  the  kinetic  equation  to  hsve  been  reported 
previously  for  exterior  flow  problems  are  those  due  to  Huang  et  ,-al . 
[7-10],  in  which  a  discrete  ordinate  method  was  applied  to  the  finite 
and  semi-infinite  flat  plate  problems . 

The  work  reported  here  is  part  of  a  coordinated  effort  by  several 
people  to  develop  numerical  procedures  capable  of  yielding  accurate 
solutions  of  the  integral  equations  for  a  sequence  of  prototype  problems 
in  kinetic  theory.  Anderson  [4,5]  has  solved  the  integral  equations  for 
several  one -dimensional  steady  state  problems,  and  Watanabe  [16,17] 
has  obtained  numerical  solutions  for  several  time  dependent  problems . 
Apart  from  their  intrinsic  interest,  these  solutions  can  be  used 
as  standards  in  assessing  the  efficacy  of  approximate  methods  for 
solving  the  Boltzmann  equation. 

In  §2,  the  integral  equations  are  formulated  in  the  context  of 
the  problem  of  a  gas  flowing  past  a  circular  cylinder.  The  free  molecule 
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form  of  these  equations  is  obtained  in  §3  and  used  in  in  the  develop¬ 
ment  of  a  numerical  method  for  solving  the  integral  equations.  Through 
the  use  of  elliptic  coordinates,  this  procedure  is  generalized  in  §5 
to  accomodate  the  case  of  an  elliptic  cylinder,  one  of  whose  axes  is 
aligned  with  the  free-stresm  direction.  It  is  also  shown  that 
singularities  develop  at  the  front  and  rear  stegnation  points  as  the 
elliptic  cylinder  degenerates  into  a  finite  flat  plate  parallel  to  the 
freestream  direction.  In  §6,  free  molecule  solutions  are  presented 
for  several  freestream  Mach  numbers  for  elliptic  cylinders  ranging  from 
a  flat  plate  parallel  to  the  freestream  direction  to  a  flat  plate 
perpendicular  to  the  freestream.  In  addition,  finite  Knudsen  number 
solutions  are  presented  for  subsonic  and  mildly  supersonic  Mach  numbers 
for  elliptic  cylinders  deviating  sufficiently  from  the  flat  plate  limit. 

2.  FORMULATION  OF  THE  INTEGRAL  EQUATIONS 

Consider  the  steady  state  flow  of  a  rarefied  gas  past  a  stationary 
circular  cylinder,  as  depicted  in  Figure  1. 

For  this  problem,  the  Krook  equation  can  be  written  in  dimensionless 
form  as  follows: 

*(vx  §  +  Vy-|)  =  ''(x){f(v;x)  -  f(v;x)l  .  (l) 

The  parameter  Tv  is  the  Knudsen  number,  v(x>)  is  a  collision  frequency 
depending  on  density  and  temperature,  and  F (v',x)  is  a  locsl  Maxwellian 
distribution 
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where  n,  q  and  T  are  density,  velocity  and  temperature,  respectively. 
These,  in  turn,  are  given  by  the  following  moments  of  the  distribution 
function: 


n(x)  =  j~  diw  f(v>;x) 

--tt  f  aVVf(V,S)  (3) 

n(x)  J 

T(x)  =  — -  f  dvtv  -  q(x)  ]2f(vjx)  . 

3n(x)  J 

In  the  foregoing  formulation,  the  following  nondimens ionalization 
conventions  have  been  implicitly  assumed.  The  characteristic  density 
and  temperature  are  taken  as  their  freestream  values  n  and  T  , 
and  the  characteristic  length  is  taken  as  the  cylinder  radius  R. 
Characteristic  values  of  velocity  and  collision  frequency  are  then 
defined  by  u  =  /kT  /m  and  v  =  v(n  ,T  ),  where  k  is  Boltzmann !s 
constant  and  m  is  the  molecular  mass .  The  Knudsen  number  is  then 
given  by  A  =  Rv/u  and  the  freestream  Mach  number  by  M  =  ''3/5  q  , 
where  q^  is  the  dimensionless  freestream  velocity. 

By  an  integration  along  its  characteristics,  the  Krook  equation 
can  be  formulated  as  a  closed  set  of  nonlinear  integral  equations 
for  the  low  order  moments  -  density,  velocity  and  temperature.  Using 
the  notation  of  Anderson  [3],  these  integral  equations  can  be  written 
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These  equations  relate  the  flaw  properties  at  a  point  x  =  (r,9)  to 
integrals  over  the  flow  field  as  viewed  from  x.  We  thus  have  a  radial 
(t)  integration  along  characteristics  passing  through  x  and  an 
angular  (<!>)  integration  over  all  such  characteristics.  For  a  given 
angle  <t>,  we  integrate  r(0  tally  along  the  characteristic 

x1  =  x  -  su  (5) 

where  u  makes  an  angle  <t>  with  the  x-axis .  Along  this  characteristic 
t  is  given  by 


(6) 


and  c  is  the  value  of  t  at  s  =  b,  the  distance  traversed  until 
a  boundary  is  intercepted.  Also,  $  is  the  angle  between  q*  and 


and  the  x-axis,  and  the  kernel  functions  are  given  by 


Kn(a,p,7)  =  Hn(a  cos  £,7)  exp  -  j 


(a  sin  ft)* 


p  <»  _ 

H  (p,q)  du  un"2  exp 

n  I2TT  J  0 


(7) 


In  the  derivation  of  equations  (4),  a  complete  accomodation  boundary 
condition  is  assumed.  Thus,  at  the  boundary  point  x  -  bu,  the  distri¬ 
bution  function  for  the  outgoing  stream  of  molecules  has  the  Maxwellian 


F(v;x  -  bu)  = 


(v  -  qj 


(2ffTB): 


where  q  and  T  are  the  velocity  and  temperature  of  the  boundary 
B  B 

at  the  point  x  -  bu  and  n  ,  a  pseudo-density  variable,  is  determined 

B 

through  the  condition  that  there  be  no  net  momentum  normal  to  the 
boundary.  In  the  context  of  the  cylinder  problem,  the  requirement  that 
the  radial  component  of  velocity  vanish  at  the  cylinder  leads  to  the 
following  integral  relation  for  n 

B 

fe  +  f 

nR(e)  = - 7=  /  a*  cos(<!>  -  9)  f  dt  n  K  (  -£•  ,<*>-  — • —  'j  ,  (9) 

^  f  ^  0  7s  n/t  / 

B  d  +  g 

the  integral  on  the  right  being  evaluated  st  the  point  on  the  surface 
to  which  9  corresponds . 

The  solid  angle  subtended  by  the  cylinder  at  the  point  (r,@) 


is  described  by  the  cone  of  influence  angles 
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Ev 


%2 


=  0  +  Arcsin 


■(10) 


For  4> €  (♦,♦_),  the  corresponding  characteristics  intercept  the 
1  &  s 

cylinder  a  distance  : 


b  =  r  cos( <t>  ••  9) 


-fi ~ 


[r  sin(<t>  -  0)  ]e 


(ID 


away.  For  4>  ,  4*^ ) ,  the  corresponding  characteristics  do  not 

intersect  the  cylinder  and  the  radial  integration  extends  over  a  s.emi- 
infinite  range.  In  general,  the  angular  integrand  is  discontinuous 

9 

> 

at  <►.  and  <t>  . 

X  d  »  , 

Equations  (4),  in  conjunction  with  (9)>  constitute  a  closed  coupled 

set  of  nonlinear  integral  equations  for  n,  q,  T,  and  n  .  In  addition, 

"  B  l 

these  integral  equations  are  singular  since  K^oc^y)1  contains'  a  term 

n-1  •,  ■ 

which  behaves  like  y  in  7  83  y  ->  0  [3].  i 

1  ! 

The  collision  frequency  v  is  commonly  chosen  as  a’  function  of 
density  and  temperature  to  match  some  continuum  transport  coefficient  ! 
of  a  particular  gas;  however,  there  is  no  unequivocal  recipe  for  this 

1 

function.  Moreover,  the  use  of  any  such  variable  collision  frequency  ' 

*  S 

model  necessitates  the  repeated  numerical  evaluation  of  the  transformation 

1  1 

(6),  which  is  embedded  in  the  radial  integration.  While  this  problem 
can  be  handled  in  a  straightforward  manner,  a  prohibitively  -large 

;  I 

computation  expense  would  be  required  for  the  class  of  problems 

1 

considered  here.  On  the  other  hand,  the  assumption  of  constant  , 
collision  frequency  reduces  the  computation  to  an  economically  tractable 

1  ! 

level.  The  emphasis  of  the  present  work  is  primarily  on  the  development 
of  numerical  methods  capable  of  yielding  accurate  solutions  of  the  kinetic 
equation,  rather  than  attempting  to  model  specific  physical  situations. 


t 

1 

I 

<5 


1 

1 

4 

'I 


1 

i 

) 

4 

% 

i 

> 

1 

! 


it 

ll 


1 


Therefore,  the  assumption  of  constant  collision  frequency  -will  be 

i  i 

.  *  !  1  : 

implicit  in  all  that  follows.  ■ 


3-  THE  FREE  MOLECULE  LIMIT  i  .  * 


The  free  polecule  solution  to  the  circular  cylinder  problem  has 
previously  been  obtained  by  Trepaud''IS:  Bruin  [13]  and  compared  with 
experimental  data.  For  cujr  purposes,  the  free > molecule  limit  furnishes 

1 

useful  information1  relevant  to  the  problem  of  solving  the  integral 
'equations  numerically.  In  the  collisionless  limit  1  =  »,  the  integral 

1 

equations  reduqe  to  the  following  closed  form  expressions : 


n  =l-‘~:  '[  d<t>  K(q  ,«,0)  +  ~  [ 

JZZ  J  *  3  00  =27r  J* 


,  3' 


/  dO  n 
27 T  B 


1  r  2  1  I  2  a 

sj  -  7=  /  id<t>  «  +  — 7=-  7  d*  “  nRs/^ 

Jo-rr  v  <t>  H  r  2iv2?r  ^  B  ^ 


3nT  t  ,nq?  =  3  +  -  “=  /  <30[K,(q  ,«,0)  +  K.(q  ,0,0)] 

v2ir  0,  ?  ^ 

!  1 


27T  J  nBTB 

*  1 


where 


1  £ 

cos‘$ 


/  ^  cos  0 


a  cos  8 
*00 


As  explicitly  recognize^  in  tl-j  above  equations,  the  free  mole  ule 

1 

solution  depends  inherently  on  the  cone  of  influence  angles  0  = 

„  b  1;2 

8  +  Arcsin(l/r ) .  Thiis  as  r  --»<»,  we  anticipate  that  the  free  molecule 
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solution  decays  to  its  freestream  limit  like  l/r,  and  that  as  r  -»1, 

the  gradients  of  the  flow  variables  approach  infinity  like  l/'/r-l. 

These  effects  stem  directly  from  the  geometric  configuration  of  the 

boundary.  As  such,  they  may  be  expected  to  occur  for  finite  Knudsen 

numbers  which  are  sufficiently  large  that  the  "boundary"  terms  (i.e. 

those  involving  n  ,  qL,  T_)  in  equations  (h)  ere  of  the  same  order 
B  B  B 

of  magnitude  as  the  integrals  over  the  interior  of  the  flow  field. 

4,  METHOD  OF  SOLVING  THE  INTEGRAL  EQUATIONS 

The  integral  equations  can  be  discretized  by  replacing  the  solution 
by  an  approximate  form  involving  finitely  rnsny  free  parameters,  replacing 
the  integrals  by  quadrature  formulae,  and  requiring  that  the  residual  - 
the  extent  to  which  the  approximate  solution  fails  to  satisfy  the  integral 
equations  -  be  minimized  in  some  sense.  For  computational  purposes,  a 
collocation  condition  was  adopted,  the  residual  being  required  to 
vanish  at  a  finite  set  of  grid  points.  The  values  of  the  flow  variables 
n,  q,  T  at  these  points  are  then  used  to  interpolate  the  flow  field 
over  the  rest  of  the  domain.  The  numerical  solution  of  the  integral 
equations  can  thus  be  considered  in  three  stages:  (l)  interpolation 
of  the  flow  field;  (2)  discretization  of  the  integrals;  and  (3)  solution 
of  the  discrete  equations. 

Since  the  flow  pattern  for  y  <  0  is  the  mirror  image  of  that 
for  y  >  0,  we  need  only  interpolate  over  the  upper  half  plane  portion 
of  the  domain;  r  >  1,  0  <  9  <  7r.  As  previously  noted,  for  sufficiently 
large  Knudsen  numbers,  the  flow  field  msy  be  expected  to  vary  roughly 
as  Arcsin(l/r)  in  the  radial  direction.  Accordingly,  the  rsdial 
coordinate  is  mapped  into  a  bounded  interval  by  means  of  the 
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transformation 


Arcs in 


(16) 


and,  as  an  additional  means  of  controlling  the  distribution  of  grid 
points,  bilinear  transformations  are  applied  to  9  and  q.  in  the 
transformed  domain  -  a  bounded  rectangle  -  the  flow  variables  sre 
interpolated  over  a  uniformly  spaced  set  of  grid  points .  A  4-point 
local  Lagrangian  scheme  is  used  in  the  radial  direction,  while  cubic 
spline  interpolation  is  used  in  the  angular  direction.  Both  inter¬ 
polation  schemes  entail  approximation  by  piecewise  cubic  polynomials; 
in  the  latter  case,  the  approximation  has  two  continuous  derivatives 
[1]. 

The  integrals  in  equations  (9)  and  (11)  are  discretized  by  means 

of  composite  low  order  Gaussian  quadrature  formulae.  Gauss -Legendre 

quadrature  is  applied  in  subintervals  in  which  the  integrand  is  analytic, 

while  special  logarithmic  and  square  root  Gaussian  formulae  are  applied 

in  subintervals  in  which  these  singularities  occur.  To  account  for 

the  discontinuities  in  the  angular  integrand  at  the  cone  of  influence 

angles  and  4>g,  separate  formulae  are  applied  to  (^,<t>2)  and 

its  complement,  each  of  these  intervals  being  allotted  a  number  of 

subinterv8ls  in  proportion  to  its  length.  A  sequence  of  arithmetically 

expanding  subintervals  is  used  for  the  radial  (t)  composite.  For 

infinitely  extensive  characteristics,  the  composite  is  applied  over 

a  finite  interval  0  <  t  <  t  .  For  t  >  t  ,  the  integrand  is 

—  —  max  max 

approximated  by  its  freestream  form,  which  can  be  integrated  directly 
using  the  relation 


'  * 
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|j;Kn(a,p,Y)  =  -K^aM)  .  (17) 

The  discretization  of  the  integral  equations  typically  yields  a 
non-sparse  system  of  several  hundred  simultaneous  nonlinear  algebraic 
equations.  These  equations  are  solved  iteratively  by  means  of  Anderson's 
Ex  rapolation  Algorithm  [2],  which  entails  a  linearization  of  the 
successive  substitution  iteration  function  about  several  of  the  most 
recent  iterates.  The  flow  variables  are  assigned  new  values  as  soon 
as  they  become  available,  as  in  a  Gauss-Seidel  iteration.  As  A 
decreases,  the  integral  equations  become  more  difficult  to  solve  numer¬ 
ically  because  the  kernel  functions  decay  more  rapidly.  Hence,  the 
best  strategy  for  obtaining  solutions  of  the  integral  equations  for 
various  values  of  A  is  to  start  with  the  free  molecule  solution  and 
proceed  in  the  direction  of  decreasing  A,  using  as  initial  iterate 
for  a  particular  value  of  A,  the  solution  corresponding  to  the  next 
larger  value. 

Another  crucial  problem  is  the  evaluation  of  the  kernel  functions. 
Since  Hn(p,q)  must  be  computed  repeatedly  during  the  course  of  solving 
the  integral  equations  numerically,  it  is  imperative  that  this  process 
be  as  efficient  as  possible.  Thus  the  computational  algorithm  is  as 
follows.  The  domain  |p|  <  2,  q  >  0  is  mapped  into  a  bounded  rectangle, 
with  a  bilinear  transformation  in  q,  and  a  rectangular  grid  of  uni¬ 
formly  spaced  points  is  laid  down  over  the  transformed  domain.  For 
n  =  2,3,4,  Hn  is  approximated  over  this  domain  by  the  bicubic  Hermite 

polynomial  which  interpolates  H  ,  its  first  partial  derivatives, 

n 

and  its  mixed  partial  derivative  st  each  of  the  grid  points .  Values 
of  Hn  for  n  >  4  are  computed  via  the  recurrence  formula 
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This  completes  our  discussion  of  the  numerical  procedure  used  to 
solve  the  integral  equations  for  the  circular  cylinder  problem.  This 
procedure  can  be  generalized  in  straightforward  fashion  to  elliptic 
cylinders,  as  described  in  the  next  section. 


5.  GENERALIZATION  TO  ELLIPTIC  CYLINDERS 

We  now  consider  a  more  general  class  of  exterior  flow  problems 
involving  flow  past  an  elliptic  cylinder  described  by 


To  facilitate  the  study  of  this  class  of  problems,  we  use  an  elliptic 
coordinate  system,  which  arises  from  the  complex  transformation  z  = 
c  cosh  p.  Equating  the  real  and  imaginary  parts  of  this  transformation 
and  making  the  appropriate  identifications,  we  obtain  the  following 
families  of  confocal  ellipses  and  hyperbolas 


In  terms  of  this  orthogonal  coordinate  system,  the  domain  exterior  to 
the  elliptic  cylinder  is  described  by 


-v  <  [x  <  ir, 


p  >  1  . 
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Equations  (9)  and  (11)  still  apply  tojthis  class  of  problems, 

:  i  ,  '  {  •  1  > 

although  the  cone  of  influence. angles  are  now  given  by  i 


*1,2  -^»*1|a4 


n>~2~ TT2^2  ~2J2 

vB  x  +  A  y  -  A  B 


(21) 


'  I  ! 

and  0  in  (11)  must  be  regarded  as  the  angle  between  the  horizontal 

■  >  !  .1  .1 
and  the  unit  outer  normal  to  the  surface.- 

I 

•  »  i  *  * 

It  can  be  shown  that  the  cone  of  the  infiuerice  of:  the  elliptic 

*  i  ,  • 

cylinder,  varies  like  l/p  as  p  -> »,  and  like  c,  +  c„  ‘Jp  -  1  as 

i  1  2  , 

p  -»1.  Thus  it  is  appropriate  to  use  the  transformation 


1}  =  Arcsin  ^ 


(22) 


|  | 

for  the .interpolation  of j the  solution.  The  rest  of  the  numerical  : 

procedure  'is  as  previously  describe^,  and  the  generalization  to  elliptic 

cylinder  is  thus  'achieved.  ■  ’  !  *  > 

'  .  ' 

It  is  instructive  to  consider ithe  behavior  of  the  free  molecule 

V 

|  I 

solution  to  the  elliptic  cylinder  problem  as  the  elliptic  cylinder  i 

,  '  '  . 
degenerates  into  a  flat  plate  parallel  to  the  freestream  direction. 

As  B  ->  0;  all  the  curvature,  of  the  elliptic  cylinder  becomes 

*  I 

concehtrated  at  the  front  and  rear  stagnation  points,  and  the  cone  of 
influence  -  and  hence  the  free  molecule  solution  -  varies  most  markedly 

v  S  v 

in  the  vicinity  of  j;hese  two  points .  In:  the  limit  B  =  0,,.  the  cone 

I  *  * 

i  ’ 

of  influence  of  the  cylinder  (in  this  case  &  flat  plate)  is  not  well-, 

1  -  1  »  *  ' 

defined  ai,  the  leading  and  trailing  edget .  Similarly,  the  limiting 

'  >  i 

'  i 

values  of  the  flow  variables  as  either  of  these  two  points  is  approached 

*  * 

1  1 
depends  upon  the  direction  of  approach. 


i  .i 


! 
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'  Free  molecule  solutions  for  a  diffusely  reflecting,  flat  plate 
have  previously  been ‘obtained  by  Ortloff  [li].  However ,  these  s'olutions 


are  tainted  by  an  incorrect  expression  for  the  distribution  function. 


The  correct  free  molecule  form  of  the  distribution  function  is 


'1 


F(v;x)  =J 


1 


(27 r) 


3/2' 


exp  - 


,V(v  -  Q Jf 


vy  0 


V. 


(27fTB) 


3^e‘ 


exp 


2T. 


B  J  ■ 


v  e  ft  , 


where  (  Q  is  the  cone  of  influence  -‘of  the  flat  plate  at  the  pc-inc  x. 
The  qualitative  behavior  of  l he  singularities  which  occur  for- 


infinite  Knudsen  number  should  also  be  manifest  at  finite  Knudsen 


■numbers,  since  within  regions  of  width  w  <<  "h  about  the  leading 
.  ,  1  ’  • 


and  trailing  edges,-  the  flow  is  essentially  1  free  molecular.  Also, 


as  B  ->  0  (and  likewise  as  A  -» 0),  the  flow  field  approaches  its 


singular  limiting  fopm,  and  the  foregoing  numerical  procedure  for 
obtaining  finite  Knudsen.'number  solutions  of  the  integral  equation 


fails  .for  "too  oblong"  elliptic  cyl.inders . 


6.  NUMERICAL  RESULTS  : 


The  integral  equations  were  solved  numerically  for  the  sets  of 


parameters  listed  in  Table  1.  Up  !to  100  grid  points  were  used  to  inter¬ 


polate  the  flow  field.  Roughly  ten  or  fewer  iterations  were  required 


-3 


to  satisfy  a  quasi-convergence  .criterion  of  10  ,  while  the  time 


taken  per  iteration  oh  the  I EM  360/65,  ranged  as  high  as1 five  minutes. 
„  .  *  *  ** 


Figures  2-5  and  6-9  depict  the  fiow  of  a  gas  past  a  circular ! 


cylinder  at  q  =1.5,,  whidh  corresponds  to  a  slightly  .supersonic 

oc  *  ,  *  1  t 


I  . 


J 


I 
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freestream  Mach  number.  The  first  set  of  contour  maps  refers  to  the 
free  molecule  limit  A  =  «,  while  the  second  refers  to  A  =  1.  Figures 
10-12  present  a  comparison  of  the  A  =  «  and  A  =  1  solutions  along 
the  rays  0-0,  v/2  and  ir.  The  similarity  of  the  two  solutions  stems 
from  the  fact  that  at  A  =  1,  the  boundary  conditions  still  have  a 
direct  influence  on  the  interior  flow  pattern.  However,  the  disturbance 
induced  by  the  cylinder  is  more  diffuse  in  the  lower  Knudsen  number 
case.  Thus,  the  square  root  singularity  at  the  cylindrical  surface 
is  less  conspicuous  at  A  =  1,  and  the  flow  variables  approach  their 
freestream  limits  more  slowly.  Also,  the  density  in  the  shielded  region 
behind  the  cylinder  is  considerably  higher  at  A  =  1. 

Figures  13-16  depict  the  flow  pattern  of  a  gas  past  an  elliptic 
cylinder  whose  semi-minor  axis  is  one-fourth  the  length  of  its  semi- 
major  axis.  Essentially  all  the  structure  of  the  flow  field  is 
concentrated  in  the  vicinity  of  the  front  and  rear  stagnation  points, 
where  the  variation  in  the  cone  of  influence  of  the  cylinder  is  greatest. 
As  the  elliptic  cylinder  degenerates  into  a  finite  flat  plate,  it 
becomes  increasingly  difficult  to  interpolate  the  solution  accurately 
in  these  two  regions . 

Figures  17-20  illustrate  the  free  molecule  solution  for  the 
aligned  flow  past  a  finite  flat  plate.  The  singularities  at  the  leading 
and  trailing  edges  are  manifested  in  the  coalescense  of  the  contour 
lines  at  these  two  points.  The  free  molecule  solutions  obtained  by 
Ortloff  [11]  do  not  exhibit  these  singularities  because  an  incorrect 
form  was  assumed  for  the  distribution  function.  Figures  21-24  illustrate 
the  free  molecule  flow  pattern  for  the  other  extreme  case  of  a  flat  plate 
perpendicular  to  the  freestream  direction.  Both  sets  of  free  molecule 
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solutions  are  for  q^  =  3,  which  corresponds  to  a  freestream  Mach 
number  slightly  greater  than  2. 

CONCLUDING  REMARKS 

.  .  ■  ..  ~  — —  ■  .  r 

The  numerical  procedure  used  to  solve  the  integral  equations  is 
most  suitable  for  high  Knudsen  number  (A  ^  l),  low  Mach  number 
(M^  <  1)  flows.  In  this  regime,  the  geometrical  configuration  of  the 
boundary  is  of  primary  importance,  and  the  flow  field  varies  essentially 
as  the  cone  of  influence  of  the  cylinder. 

At  lower  Knudsen  numbers,  the  "width"  of  the  kernel  functions 
decreases,  and  the  effect  of  the  boundary  conditions  on  the  interior 
flow  pattern  is  less  direct.  More  precisely,  the  boundary  integrals 
in  (4)  at  a  fixed  point  x  decrease  in  magnitude  as  \  -»0.  Thus  a 
numerical  method  tailored  to  the  variation  of  the  cone  of  influence 
angles  is  not  necessarily  appropriate  for  solving  the  integral  equations 
for  low  Knudsen  numbers. 

At  higher  Mach  numbers,  the  gradients  in  the  flow  variables  become 
larger,  posing  additional  numerical  difficulties.  Due  to  the  shielding 
effect  of  the  cylinder,  density  may  vary  by  several  orders  of  magnitude. 
For  example,  the  free  molecule  density  at  the  rear  stagnation  point 
is  of  order  l/q^  exp  -(q^/2)  as  q  -> «  (this  result  can  be  obtained 
from  equations  (l4)  and  (15)).  Also,  at  finite  Knudsen  numbers,  the 
occurrence  of  shock  waves  presents  an  added  complication.  The  question 
of  how  to  handle  these  effects  numerically  remains  unresolved. 

Also  unsettled  is  the  question  of  how  to  adapt  the  integral  equation 
approach  to  the  flat  plate  problem.  It  would  seem  that  the  singularities 
at  the  leading  and  trailing  edges  should  be  taken  account  of  numerically, 


yet  it  -is  unclear  hew  this  should  be  accomplished.  At  any  rate,  one 
cannot  assume  a  priori  that  the  effect  of  these  singularities  is 
negligible. 
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FIG.  10  DENSITY  PROFILES 
Ocn  =  1.5.  T„  =  l 
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FLOW  PAST  AN  ELLIPTIC  CYLINDER 
X*l  0^=1  Tb=1 
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DW  PAST  A  FINITE  FLAT  PLATE 
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